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We investigate the possibility of forming deeply bound ultracold RbCs molecules by a two-color 
photoassociation experiment. We compare the results with those for Rb2 in order to understand the 
characteristic differences between heteronuclear and homonuclear molecules. The major differences 
arise from the different long-range potential for excited states. Ultracold *^Rb and ^^^Cs atoms 
colliding on the X^E+ potential curve are initially photoassociated to form excited RbCs molecules 
in the region below the Rb(5S) 4- Cs(6Pi/2) asymptote. We explore the nature of the = 0^ 
levels in this region, which have mixed A^E"*" and b'^fl character. We then study the quantum 
dynamics of RbCs by a time-dependent wavepacket (TDWP) approach. A wavepacket is formed by 
exciting a few vibronic levels and is allowed to propagate on the coupled electronic potential energy 
curves. We calculate the time-dependence of the overlap between the wavepacket and ground-state 
vibrational levels. For a detuning of 7.5 cm~^ from the atomic line, the wavepacket for RbCs 
reaches the short-range region in about 13 ps, which is significantly faster than for the homonuclear 
Rb2 system; this is mostly because of the absence of an _R~^ long-range tail in the excited-state 
potential curves for heteronuclear systems. We give a simple semiclassical formula that relates the 
time taken to the long-range potential parameters. For RbCs, in contrast to Rb2, the excited- 
state wavepacket shows a substantial peak in singlet density near the inner turning point, and this 
produces a significant probability of deexcitation to form ground-state molecules bound by up to 
1500 cm~^. The short-range peak depends strongly on nonadiabatic coupling and is reduced if 
the strength of the spin-orbit coupling is increased. Our analysis of the role of spin-orbit coupling 
concerns the character of the mixed states in general and is important for both photoassociation 
and stimulated Raman deexcitation. 



I. INTRODUCTION 

The study of cold molecules, below 1 K, and ultracold 
molecules, below 1 mK, offers many new opportunities 
for chemical and molecular physics. Cold molecules open 
up a new regime for molecular collisions in which clas- 
sical physics breaks down completely and a quantal de- 
scription is needed for all degrees of freedom. Collisions 
in this regime are dominated by long-range forces and 
exhibit resonance phenomena that can be controlled by 
applied electric and magnetic fields. Cold molecules also 
present new opportunities for precision measurement and 
open up possibilities for measuring quantities such as the 
dipole moment of the electron and the time-dependence 
of fundamental "constants" . 

There is particular interest in ultracold polar 
molecules, because dipolar species interact more strongly 
and at much longer range than non-polar species. Dipo- 
lar quantum gases are predicted to exhibit new phenom- 
ena such as anisotropic Bosc-Einstcin condensation [l| 
and may have applications in quantum information pro- 
cessing [3 . They also provide opportunities for engineer- 
ing highly correlated quantum phases Q . 

Molecules can be formed in ultracold atomic gases by 
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both magnetoassociation and photoassociation [1, 0, @1 , 
and there have been considerable advances in using such 
methods to produce heteronuclear alkali metal dimers 
0, i, i, E El, [13, m, [ll- However, both photo- and 
magnetoassociation initially produce molecules in very 
high vibrational levels, which are only very weakly dipo- 
lar even for heteronuclear species. There is therefore 
great interest in producing ultracold molecules in low- 
lying vibrational states. For example. Sage et al. [1, Q 
have succeeded in producing small numbers of ultracold 
RbCs molecules in the ground vibronic state by a four- 
photon photoassociation scheme using continuous-wave 
lasers, while Deiglmayr et al. have produced LiCs 
molecules with a two-photon scheme. Both these ap- 
proaches include a spontaneous emission step, but Ni et 
al. |l3 have very recently produced KRb molecules in the 
lowest vibrational levels of the lowest singlet and triplet 
electronic states by magnetoassociation followed by a co- 
herent two-photon process (stimulated Raman adiabatic 
passage, STIRAP). Danzl et al. [1^1 have carried out 
analogous experiments on the homonuclear molecule Cs2 , 
while Lang et al. [13] have produced Rb2 molecules in 
the lowest level of the lowest triplet state. 

Heteronuclear molecules differ from homonuclear 
molecules in several ways. The most important differ- 
ence is that the excited-state potential curves correlating 
with -I- ^P atoms have an behavior at long range 
in the heteronuclear case but an behavior in the 
homonuclear case, due to the resonant dipole interaction. 
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FIG. 1: The potential curves and laser pulse scheme consid- 
ered here. The potential curves correspond to the ground 
X^E^ and the excited A'^S"'" and b"^n states. A picosecond 
laser pulse is used to excite few vibronic levels in the excited 
state to form a wavepacket. A second laser pulse is used 
to dump the wavepacket to the ground X^E"*" state after an 
appropriate time delay. The vibrational wavefunction of the 
target level of the dump pulse, v = 54, is also displayed. 



Because of this, the Franck-Condon factors for photoas- 
sociation are quite different. In addition, the density of 
vibrational levels in the electronically excited state is dif- 
ferent for and potentials. 

Magnetoassociation must be carried out in tight traps 
at very low temperatures. However, there would be ad- 
vantages in producing molecules in low-lying states at 
the somewhat higher temperatures (in the microkclvin 
regime) that are available in magnetooptical traps 
(MOTs), where the number of atoms available is of- 
ten much larger. A very promising approach for this is 
photoassociation using shaped laser pulses. Short-lived 
molecules are formed in the excited electronic state by 
photoassociation using a short laser pulse (pump pulse) 
during the collision of two ultracold atoms. The excited 
molecules are then stabilized by stimulated emission us- 
ing a second laser pulse (dump pulse) into the bound 
vibrational levels of the ground electronic state. Luc- 
Koenig et al. [H, ^ and Koch et al. |20l . l2l| have sim- 
ulated this process for homonuclear diatomic molecules 
such as Cs2 and Rb2, and initial experiments on these 
systems have been reported by Salzmann et al. [13, [23| 
and Brown et al. (2^ . Similar experiments can be envi- 
sioned for heteronuclear molecules, with the use of evo- 
lutionary algorithms or other strategies to maximize the 
rate of production of ground-state molecules. However, a 
theoretical study of short-pulse photoassociation of het- 
eronuclear molecules has not yet been carried out. In 
this paper, we consider RbCs as a prototype heteronu- 
clear molecule and we study the possibility of forming 
ground-state RbCs molecules through pump-dump pho- 
toassociation. 

We consider the photoassociation of a colliding pair of 
®^Rb and ^'^^Cs atoms, initially in their ground state. 
Absorption of a photon red-detuned from the atomic 



Rb(5S) -I- Cs(6Pi/2) asymptote forms an excited-state 
molecule around internuclear distance Rl, as shown in 
Fig. [TJ When spin-orbit coupling is included there are 8 
electronic states that correlate with the excited 5Si/2 + 
6P1/2 and 5S1/2 + 6P3/2 asymptotes: two 0^ states, two 
0^ states, three $7 = 1 states and one il = 2 state [25| . 
These correlate at short range with the A ^2+, b "^H, B ^H 
and c'^S"'" electronic states in Hund's case (a) labelling. 

Absorption of a photon from the ground X^I]+ (0+) 
state can produce excited-state molecules in O"*" or 1 
states, while absorption from the a3i;+ state (with 0" 
and 1 components) can produce 0+, 0", 1 and 2 states. 
We focus here on the 0"'' excited states, which may be 
formed by combining the A^S+ and b'^H states. The 
c'^E+ state does not contribute because it has only 0~ 
and 1 components and the B ^H state does not contribute 
because it has only an Q = 1 component. Rotational 
couplings that would connect states of different Q are ne- 
glected. This 2-state approximation is analogous to the 
approach used for the 0+ states in Cs2 and Rb2 [1^ , for 
which the strongly mixed singlet-triplet character of the 
vibronic wavefunctions was found to produce enhanced 
formation of ground-state molecules with binding ener- 
gies on the order of 10 cm~ 

ana- 

The pump pulse produces a non-stationary state made 
up of long-range excited-state levels of mixed singlet and 
triplet character. The corresponding wavepacket prop- 
agates towards short range under the influence of the 
excited-state potentials. After a suitable time delay, 
when a sufficient amount of the wavepacket has reached 
the short-range region, a dump pulse is activated to 
transfer the molecules into vibrational levels of the X ^S]+ 
electronic ground state. 

We simulate the entire pump-dump process using a 
time-dependent wavepacket (TDWP) approach. Our 
goal is to optimize the parameters of the pump and dump 
pulses to maximize the production of molecules in deeply- 
bound levels of the ground state. We find that a crucial 
factor is the magnitude of the spin-orbit coupling near 
the avoided crossing. If the spin-orbit coupling is suffi- 
ciently strong, then the levels that lie on the lower adi- 
abatic curve, which correlates with 5S + 6P1/2 at long 
range, have mainly triplet character at short range. Such 
levels have low intensities for deexcitation. We therefore 
form the excited-state wavepacket by selective photoas- 
sociation into vibronic levels of the A ^E+- b '^H electron- 
ically excited states that are strongly mixed by nonadi- 
abatic coupling and have enhanced singlet character at 
short range. We estimate the parameters required for the 
dump pulse and the time delay by analyzing the time- 
dependence of the overlap between the singlet part of 
the wavepacket and the vibrational levels of the ground 
electronic state. 

We find an important difference between photoassoci- 
ation for heteronuclear and homonuclear molecules. The 
excited states of heteronuclear molecules have a lower 
density of vibrational levels near dissociation than those 
of homonuclear molecules, because of the long-range 
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term in the excited-state potentials for homonu- 
clear species. Because of this, the atoms of heteronuclear 
molecules produced at long range experience a signifi- 
cantly faster classical acceleration towards one another 
than homonuclear molecules created at the same bind- 
ing energy. The time delay required between the pump 
and dump pulses is thus smaller in the heteronuclear case 
(about 13 ps for RbCs with a detuning of 7.5 cm"-'^ from 
the atomic line) than in the homonuclear case. However, 
even for heteronuclear molecules, there is little probabil- 
ity of producing ground-state molecules within the dura- 
tion of a single pulse of less than a few ps. 

The paper is organized as follows. In section II, we dis- 
cuss the potential energy curves and spin-orbit coupling 
functions used in this study. In sections III and IV we 
explore the dependence of the zero-energy s-wave scatter- 
ing lengths on the interaction potentials and discuss their 
influence on the Franck-Condon factors for photoassocia- 
tion. In section V we describe wavepacket studies of the 
photoassociation process and the subsequent formation 
of ground-state molecules. At each stage, the results for 
RbCs are compared with those for Rb2 in order to estab- 
lish the similarities and differences between heteronuclear 
and homonuclear molecules. 
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FIG. 2: The diagonal and off-diagonal components VF™(ii) 
and Wi^{R) of the spin-orbit coupling as a function of R. 
The avoided crossing between the two adiabatic potentials is 
shown in the inset. 



for the A and b states respectively. 

Spin-orbit coupling between the A and b ■^11 states 
plays an important role. Following Marinescu and Dal- 
garno [s^] and Aubert-Frecon et al. [3^ , the effect of the 
spin-orbit coupling for the 0+ states is included in terms 
of a 2x2 electronic Hamiltonian matrix. 



II. POTENTIAL CURVES 

Several studies of the potential energy curves for RbCs 
have been published. Obtaining accurate potentials is 
not a straightforward task, because electronic structure 
calculations are very difficult for such heavy atoms. Ex- 
perimental results are available for some states [28l . [29j , 
but are sparse or nonexistent for others. Only some of 
the theoretical potentials include spin-orbit coupling ef- 
fects; for example the relativistic curves of Kotochigova 
and Tiesinga [2g include avoided crossings which are ab- 
sent in the older pseudopotential results of AUouche et 
al. [13 . 

The potential curves for the electronic states of 
RbCs considered here, neglecting spin-orbit coupling, are 
shown in Fig. [TJ From the viewpoint of photoassociation, 
the most important parts of the curves are the long-range 
tails. In this study we have chosen for simplicity to use 
the ah initio results of AUouche et al. [soj at short range 
(bond lengths R < 15 ag). For the long-range potentials, 
in the ground state, we use the experimentally derived 
parameters of Fellows et al. [2^ , with Ce replaced by the 
more recent value of Derevianko et al. [3l|. These pa- 
rameters are equivalent to column V of Table 2 in Rcf. 
[3^ . The Lc Roy radius for the ground X ^S]+ state is 
24.4 oq, but a smooth match between the long-range and 
short-range potentials was best achieved at 16.5 ag. For 
the excited singlet and triplet states, we use the theoreti- 
cal long-range parameters of Marinescu and Sadeghpour 
[33| . corresponding to the Rb(5S) -|- Cs(6P) asymptote. 
For the excited states the matching between long-range 
and short-range curves was achieved at 16.5 ao and 24 qq 
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W^^iR), 
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where VAis+(i?) and Vb3n(-R) are the Born- 
Oppenheimer potentials neglecting spin-orbit coupling. 
The diagonal and off-diagonal spin-orbit coupling ele- 
ments are denoted W™{R) and W^q{R) respectively 
and are functions of the internuclear separation R as 
shown in Fig. [21 At long range. 
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where AE'qq is the spin-orbit splitting of the P state of 
atomic Cs. Bergeman et al. |36| give a non-zero matrix 
element connecting the = components of the h^H 
and c'^S"'" states but this is appropriate only for the 0~ 
component and does not influence the 0+ states. We use 
the spin-orbit coupling functions obtained by Fellows and 
Bergeman (37| . which were obtained by fitting to results 
from Fourier-transform spectroscopy [231 . The diagonal 
and off-diagonal spin-orbit couplings show a significant 
dip from their asymptotic values near the crossing be- 
tween the A^S''' and b'^II potential energy curves. The 
two adiabatic potentials obtained by diagonalizing the 
electronic Hamiltonian of Eq. ([T]) are shown as an insert 
in Fig. [1 

For Rb2 the potentials are more accurately known. In 
the present work, the ground-state X^E+ potential is 
obtained by combining the short-range results of Seto et 
al. [3^ with the long-range coefhcients of Marte et al. 
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FIG. 3: RbCs s-wave scattering length for the X state, as 
a function of a scaling factor used to multiply the potential. 



[39|. The excited-state potentials arc taken from Berge- 
man et al. [40|. who obtained potential curves and spin- 
orbit coupling functions for the A^I]+ and b'^IIu states 
by combining photoassociation spectroscopy with short- 
range ab initio results [4l| . The corresponding curves for 
Rb2 (not shown here) are qualitatively similar to those 
for RbCs, except for the important difference in long- 
range behavior. 



III. SCATTERING LENGTH 
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FIG. 4: Upper panel: Average bond length (R) for vibrational 
levels of RbCs in its A ^E"*" state as a function of the binding 
energy. Lower panel: Franck-Condon factors for photoassoci- 
ation (in arbitrary units), calculated for a ground-state scat- 
tering length a — 66.3 ao- Note that the largest factors are 
off the top of the scale in this figure. 



A crucial parameter affecting low-energy scattering 
properties and the positions of high- lying excited states is 
the s-wave scattering length a. Calculations by Jamieson 
et al. [sl] produced values for a ranging from 380 to 1 oq 
for the X state of RbCs, depending on the subtleties 
of the potential parameters. There are further uncertain- 
ties due to the choice of hyperfine states involved. In 
addition, even if the potential is known accurately, there 
is the possibility of controlling the scattering length with 
applied fields. We therefore investigate the dependence 
of photoassociation on the scattering lengths for both the 
ground and excited states. 

Scattering lengths were calculated by solving the 
Schrodinger equation for zero-energy s-wave scattering 
numerically using Numerov integration. The long-range 
wavefunction was matched to Bessel and Neumann func- 
tions of fractional order [l^l to take account of the long- 
range potential and avoid the need to propagate to 
excessively large distances. We adjusted the scattering 
length by scaling the whole potential curve by a con- 
stant factor. As is well known, the scattering length is 
extremely sensitive to such scaling and passes through a 
pole whenever there is a bound state at exactly zero en- 
ergy. For ^■^Rb^^'^Cs a complete cycle is achieved within 
a scaling of just over 1%, as can been seen from Fig. [31 



Rb2 shows similar behavior. 



IV. FRANCK-CONDON FACTORS AND THEIR 
DEPENDENCE ON SCATTERING LENGTHS 

We first calculate the Franck-Condon factors (FCFs) of 
relevance to photoassociation experiments for ^^Rb^^'^Cs 
and compare the results with those for the homonuclear 
^^Rb2 system. The purpose of this Section is to un- 
derstand the qualitative differences between RbCs and 
Rb2, which stem from the different long-range tails of 
the excited-state potential (proportional to for RbCs 
and for Rb2). These qualitative differences are most 
simply illustrated using calculations that neglect the cou- 
pling between the A and b'^II states, and that is the 
approach we use in this Section. However, coupled cal- 
culations of FCFs have been carried out by Luc-Koenig 
[isj , and indeed the calculations described in Section lV CI 
to IV El below use fully coupled wavefunctions. 

The FCFs for photoassociation are calculated as the 
squares of overlap integrals between bound vibrational 
states in the excited A^I]+ electronic potential and a 
low-energy scattering wavefunction on the ground-state 
potential. Both the bound and scattering states were 
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FIG. 5: Franck-Condon factors for RbCs photoassociation, 
calculated for a variety of ground-state scattering lengths. 
The scale is arbitrary but is consistent between graphs. 



obtained by Coolcy-Numerov propagation [4^ for rota- 
tional angular momentum = 0. A LeRoy-Bernstein 
analysis [45| confirmed that every bound state had been 
found. The s-wave scattering function was calculated at 
a collision energy of 1 mK and normalized by setting its 
maximum amplitude to 1 as i? ^ oo. The ground-state 
potential docs not die off to 1 mK until i? '--^ 110 ao, so 
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FIG. 6: Upper panel: Average bond length (R) for vibrational 
levels of Rb2 in its A state as a function of the binding 
energy. Lower panel: Franck-Condon factors for photoassoci- 
ation. The scale is arbitrary but is consistent with that for 
RbCs in Figs. |4] and [S] Note that the largest factors are off 
the top of the scale in this figure. 



that the relative values of the resulting Franck-Condon 
factors should be valid for excitation to states dominated 
by distances R < 100 qq. 

Figure [3] shows the calculated FCFs for RbCs with 
a = -1-66.3 oq as a function of binding energy (detuning 
from the Cs atomic line), and the corresponding values 
of (R) for the excited vibrational states. The FCFs show 
a typical oscillating pattern as the bound and continuum 
states come into and out of phase with one another. For 
RbCs each peak of the envelope encompasses around 4 
vibrational levels of the excited state. The most intense 
peak is just below the dissociation threshold, and the 
state at the center of this peak has (i?) sa 71 qq. The 
next peaks in order of decreasing intensity correspond to 
(R) w 47, 37, 32, 29 and 26 aQ. The largest FCF occurs 
for the least-bound state, with (R) = 150 ao, which is 
not part of an envelope. However, this state is bound by 
only 0.0007 cm~^, and it is likely that frequencies this 
close to the atomic line would be blocked in an actual 
experiment. 

To investigate the effect of adjusting the ground-state 
scattering length, we have repeated the Franck-Condon 
calculations for a = 135.8, 33.5, and —68.8 ao, corre- 
sponding to potential scaling factors of 0.9965, 1.0015 
and 1.0039 respectively. The results are shown in Fig. 
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[5j It may be seen that changing the scattering length 
shifts the positions of the peaks in the envelope of the 
intensity distribution. Nevertheless, the peak intensi- 
ties themselves follow a curve that is only very weakly 
potential-dependent. The peak intensities are thus al- 
most a single- valued function of binding energy. 

Scaling the potential curve for the excited state, while 
keeping the ground-state scattering length constant, 
shifts the positions of the individual lines but does not 
alter the positions of the peaks in the intensity distri- 
bution. This effect is shown for a — -1-66.3 Uq in Fig. [8l 
the frequencies that correspond to vibrational levels shift 
as the potential is scaled, but the envelope of the FCFs 
does not change. Since a laser pulse will always cover a 
range of binding energies, it is this envelope that matters 
more than the specific position of the levels and the over- 
all transition probabilities will not be greatly affected by 
changes in the excited-state scattering length. 

The Franck-Condon factors for photoassociation to 
form Rb2 are shown in Fig. [6l The potential curve used 
here has a = 157.9 oq. Once again the FCFs show an 
oscillatory structure. The major difference from RbCs 
is that the oscillations in the envelope of line intensities 
are slower as a function of binding energy and the lines 
themselves are much more densely packed because of the 
greater density of states for an potential. The over- 
all values of the FCFs are also somewhat larger in the 
homonuclear case. 

Another significant difference between the heteronu- 
clear and homonuclear cases is shown in Fig. [71 On a 
log-log plot, the intensities of the peak envelopes vary 
almost linearly with binding energy Ed — near disso- 
ciation. This corresponds to a power-law dependence on 
binding energy. However, the powers involved are clearly 
different for the two cases: about [E-q — Ey)~^ for the 
homonuclear case and {ED—E^)~^^^ for the heteronuclear 
case. This will be significant in designing pulsed-laser 
experiments, because the intensity will fall off consider- 
ably faster with detuning (from the atomic line) in the 
homonuclear case than in the heteronuclear case. Thus a 
broad pulse will excite a wavepacket in which the relative 
population of deeper levels is lower in the heteronuclear 
case. 



V. TIME-DEPENDENT CALCULATIONS 



A. Vibrational periods 



In section IV Bl we study the photoassociation process 
in a time-dependent wavepacket approach. However, we 
can make a rough estimate of the time needed for the 
wavepacket to evolve from long to short range by consid- 
ering the vibrational half-period. 
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FIG. 7: Franck-Condon factors for different values of the 
ground-state scattering length, shown on a log-log plot for 
RbCs (upper panel) and Rb2 (lower panel). Note that the 
peak values lie on a straight line in each case, but with a 
different slope, corresponding to a different power-law depen- 
dence on binding energy. 



where v is the classical vibrational frequency. The vibra- 
tional spacing dEy/dv may be calculated exactly for a 
particular potential curve but it is instructive to consider 
its value in long-range theory. If only the leading term in 
the long-range potential is included, V{R) = — C„i?~", 
the vibrational period may be written semiclassically [46j 



^^K,,{ED-Ey)(-+'y' 
dv 



where 



(27r)i/2/i / nT{l + 1/n) 



y.^l^&J'' ^i/2cy" Vr(l/2 + l/n) 



(3) 



(4) 



/X is the reduced mass and F is the gamma function. Le 
Roy has tabulated numerical values of the coefficients 
pq . The times obtained for wavepackets corresponding 
to each of the main peaks in Fig. 2] are given in Table |T] 
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FIG. 8: Franck-Condon factors for different scalings of the up- 
per potential with the unsealed ground-state potential. Black, 
blue and red symbols are for scalings of 1.0, 1.001 and 1.005 
respectively. 



TABLE L Vibrational half-periods for an RbCs molecule at 
energies corresponding to maxima in the intensity distribu- 
tion for photoassociation. The index m labels the successive 
maxima, starting at the dissociation limit. 
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vary from 


over 1 



at small detunings to under 1 ps for a detuning of —315 
cm~^. Although the positions of the individual maxima 
are sensitive to details of the potential, the time taken 
is an almost single- valued function of the binding energy 
(laser detuning) for a particular long-range potential and 
reduced mass. 

The times taken for RbCs are significantly shorter than 
those for Rb2, shown in Table HIl For example, for a 
detuning of 25 cm^^ from the atomic line, the classical 
time is about 13 ps for Rb2 but 5 ps for RbCs. 

TABLE II: Vibrational half-periods for an Rb2 molecule at 
energies corresponding to maxima in the intensity distribu- 
tion for photoassociation. The index m labels the successive 
maxima, starting at the dissociation limit. 



m 


binding energy 


(7?) 


relative 


ii/2 




(cm-i) 


(«o) 


intensity 


(ps) 


1 


-0.017 


500 


12 


6300 


2 


-0.038 


380 


8.5 


3200 


3 


-0.10 


270 


5.2 


1400 


4 


-0.49 


160 


2.4 


380 


5 


-5.1 


75 


0.36 


54 


6 


-27 


43 


0.043 


13 


7 


-56 


34 


0.014 


7.3 


8 


-91 


29 


0.007 


4.9 


9 


-135 


26 


0.004 


3.5 



The classical model used here is of course considerably 
oversimplified and neglects features such as the spreading 
of a wavepacket as it propagates inwards. A full treat- 
ment of the time evolution requires quantum-mechanical 
calculations as described below. 



B. Wavepacket model and methods 

In the time-dependent wavepacket approach, the pho- 
toassociation reaction starts with two cold atoms in 
the ground electronic state with relative kinetic energy 
i?kin- Absorption of a photon red-detuned by an energy 
fkOL from the atomic resonance line at hcunt produces a 
molecule in the excited electronic state as shown in Fig. 
[T] The Hamiltonian describing optical transitions be- 
tween the X ^E"'" electronic ground state and the coupled 
A^E+ and b^II excited states can be represented in the 
diabatic basis as 



H 



ti{R) ■ E{t) t 
V V2Wi^{R) 



fi{R) ■ E{t) 

Vp^iY,+ {R) - fy^L 



\ 

V2Wi^{R) 
V^m{R)-Wi^{R)-nwL I 



(5) 



where T is the kinetic energy operator and Vi{R) are the and rotating- wave approximations, the coupling between 
respective potential energy curves. Assuming the dipole 
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the X-^E+ and A-'^S+ electronic states is represented by 
the scalar product between the transition dipole moment, 
and the polarization vector of the laser field, E{t). 
In the present work the i?-dependence of is ne- 

glected and it is represented by its asymptotic value. 
The A^E+ and b'^II excited states are coupled by the 
spin-orbit interaction. The hyperfine interaction can be 
neglected because the associated timescale is much longer 
than the femtosecond or picosecond laser pulses, so that 
the hyperfine interaction is not resolved for the processes 
considered in the present study. 

The Hamiltonian, Eq. ([5]), is represented on a Fourier 
grid with a variable grid step [l^ . Details of the mapped 
Fourier Grid Hamiltonian (FGH) method are described 
in refs.H^ and l48l. Using the FGH method, we are able 
to extend the spatial grid to i? 10'^ gq using only 1023 
grid points; this is sufficient both to represent the bound 
vibronic levels and to approximate the scattering contin- 
uum in the ultracold temperature regime by box states 

M- 

To describe the short-pulse photoassociation pro- 
cess, we solve the time-dependent Schrodinger equation 
(TDSE), 

ih-^^^'it) = H(t)*(i), (6) 

by expanding the evolution operator exp[—ifit/h] in 
terms of Chebychev polynomials [s^l. Because of the 
large extent of the grid and the very small kinetic energy 
of the system, it is not necessary to enforce an absorbing 
boundary condition at the edge of the grid. 



C. Analysis of the coupled vibronic wavefunctions 

Diagonalization of the Hamiltonian of Eq. ([5]) with 
E{t) set to zero yields the binding energies and vibronic 
wavefunctions and allows the calculation of FCFs, rota- 
tional constants, etc. For RbCs it produces 136 bound 
vibronic levels for the X^E^ ground electronic state and 
398 levels for the coupled A^E+ and b^H electronic 
states. The eigenfunctions of the excited electronic states 
are perturbed by the spin-orbit coupling and have mixed 
singlet and triplet character. The highest excited-state 
level has a binding energy of 0.002 cm~^, which differs 
slightly from the result in Section ITvl because the latter 
neglected the coupling between the singlet and triplet 
states. 

When two electronic states interact and their curves 
cross, there are two quite different limiting cases that may 
be considered to be "uncoupled" . If the coupling near the 
crossing point is very weak (compared to the local vibra- 
tional spacings), the crossing is only weakly avoided and 
the eigenstates of the coupled system are close to those of 
the individual diabatic (crossing) electronic states. Con- 
versely, if the coupling is strong, the crossing is strongly 
avoided and the best zeroth-order picture is to consider 



eigenfunctions of the adiabatic states, defined by the non- 
crossing upper and lower adiabatic curves. In this case 
the adiabatic states themselves change character from 
one side of the crossing to the other, so that in the present 
case the zeroth-order states are predominantly singlet on 
one side of the crossing and predominantly triplet on the 
other. The couplings between the adiabatic states are 
provided by nonadiabatic couplings (which are related to 
off-diagonal matrix elements of d/dR and d'^ /dB?). 

The wavefunction for a mixed vibronic state may be 
written 

^n^)V^l+V3(^)^'3 (7) 

where the functions ^\{R) and ^^{R) are vibrational 
wavefunctions and tj}! a-nd -03 are electronic functions for 
the singlet and triplet excited states respectively. It is 
useful to compare the eigenstates that are obtained for 
two different cases, as shown in Fig. [51 The two panels on 
the left show RbCs eigenstates calculated with the spin- 
orbit coupling functions fixed at their asymptotic value, 
^^so = 554.04 cm-i (corresponding to W^^ = 184.68 
cm~^, which is considerably larger than the experimen- 
tally derived value [sj around the crossing point of the 
diabatic curves {W^q « 100 cm^^, see Fig. [S]). These 
will be referred to as ASO (asymptotic spin-orbit) calcu- 
lations. The two levels shown are for vibrational indices 
v' = 359, with binding energy 22.19 cm~^, and v' = 365, 
with binding energy 10.01 cm~^. It may be seen that 
for both ASO wavefunctions there is a switchover from 
singlet to triplet character around the avoided crossing 
[R ^ 10 ao), but in opposite directions. This occurs be- 
cause the states are nearly adiabatic: the function for 
v' = 365 lies on the upper adiabatic curve, so has mainly 
singlet character inside the crossing and mainly triplet 
character outside. The function for v' ~ 359 lies on 
the lower adiabatic curve, so that it has mostly triplet 
character at short range. Outside the crossing, it has 
mainly singlet character in the region where the singlet 
and triplet curves are separated by more than the spin- 
orbit coupling. At long range (i? > ~ 20 ao), however, 
it reacquires the triplet admixture characteristic of an 

= state at the 5S + 6P1/2 threshold, yielding a 2:1 
ratio of triplet to singlet probability. It stretches to con- 
siderably greater internuclear distance because of the dif- 
ferent turning points for the two adiabatic states. 

The near-adiabatic vibronic levels obtained for the 
ASO case provide a context to understand the eigenstates 
obtained with the i?-dependent spin-orbit (RSO) cou- 
pling function, shown in the panels on the right of Fig.[9l 
In the RSO case the levels with index numbers v' = 359 
and 365 have binding energies 17.79 cm~^ and 7.53 cm~^, 
respectively. The v' = 359 level exists mostly on the lower 
adiabatic curve, and has triplet character with only a 
small singlet admixture at short range and switches over 
to singlet character (with the same characteristic triplet 
admixture as before) outside R = 10 oq. Nevertheless, 
there is significant nonadiabatic mixing which introduces 
upper-state character, shown by the peaks in probabil- 




FIG. 9: Eigenfunctions of the coupled excited-state potentials corresponding to u = 359 and 365. Left-hand panels: results 
obtained with constant asymptotic spin-orbit coupling (ASO). Right-hand panels: results obtained with 7i-dependent spin-orbit 
coupling (RSO). 



ity densities around the outer turning point of the upper 
curve near i? ~ 17 Oq. The v' = 365 level is more strongly 
mixed, with a larger contribution from the upper adia- 
batic state that makes it predominantly singlet at short 
range. However, this state too has substantial lower- 
state character so that it does have density out to the 
outermost turning point {R ~ 33 ao in this case). It is 
notable that both RSO eigenstates have significantly en- 
hanced singlet density near the outer turning point of the 
upper curve (R ~ 17 ao). This is the feature that in Rb2 
was responsible for enhanced deexcitation to levels of the 
ground state bound by up to 10 cm~^ Ill|,[13l- However, 
the wavefunction for v' = 365 also has substantial singlet 
density near the inner turning point of the singlet state; 
this is a feature that was not observed for Rb2 (see Fig. 
2 of ref. nil), and as will be seen below it allows deexci- 
tation to even more deeply bound levels. The fact that 
the wavefunction can have significant singlet density both 
at this turning point and at the outermost turning point 
arises from nonadiabatic coupling between the upper and 
lower adiabatic states. 



In an alternative approach, the nonadiabatic mixings 
could be shown in an adiabatic representation, display- 
ing the components of the vibrational wavefunctions on 
the two states. This would emphasize the nearly adia- 
batic character of the ASO wavefunctions but would hide 
their switchover from singlet to triplet character near the 
avoided crossing. 

The character of the mixed levels, as a function of vi- 
brational quantum number, may be seen very clearly in 
the calculated rotational constants for the coupled lev- 
els (By ~ {h/{8n'^cfiR'^))), which arc shown in Fig. [10] 
for both the ASO and RSO cases. In the ASO case, the 
rotational constants form two almost independent series, 
with the higher values corresponding to levels that are 
predominantly on the upper adiabatic curve. The up- 
per state has a larger vibrational spacing than the lower 
state. In the RSO case, the peaks in rotational constants 
correspond to levels with a significant contribution from 
the upper adiabatic state. However, it may be seen that 
in this case the nonadiabatic coupling is strong enough to 
spread the character of each upper-curve vibronic state 
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FIG. 10: Rotational constants B(, — {h/ {?,-k^ c^B?)) corre- 
sponding to the vibronic levels of the coupled A^E+ + 6^11 
excited electronic states for RbCs. 
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FIG. 11: Rotational constants B'^ — {h/ {^-k"^ c^B?)) corre- 
sponding to the vibronic levels of the coupled A^E^ + h'^Il 
excited electronic states of Rb2. 



across several eigenfunctions of the coupled problem, so 
that the peaks are not very large. There are no eigen- 
states of the coupled problem that are predominantly on 
the upper adiabatic curve. 

The rotational constants for RbCs may be compared 
with those for Rba [H, SO], shown in Fig. [Ill Rba ex- 
hibits nonadiabatic mixing that is intermediate between 
the ASO and RSO cases for RbCs. For Rb2 the extent of 
the mixing has been shown to be isotope-dependent [5l| . 

The existence of bound states with singlet character 
both at long range and near the inner turning point will 
also be important in designing fixed-frequency schemes to 
populate deeply-bound ground-state levels by STIRAP 
and related methods. Such schemes will be most efficient 
when they proceed via mixed levels with population on 
both the upper and lower adiabatic curves. 



TABLE III: The parameters for the pump and dump pulse: 
detuning Sf^ from the atomic line, intensity II, temporal 
width tl, spectral bandwidth 5^^, and integrated pulse energy 
per area. 



pulse if 


II 




(5^/(27rc) 


Energy /Area 


(cm-') 


(W cm" 


(ps) 


(cm-i) 


(J m-^) 


pump (ASO) -10.01 


16.86 


5.00 


2.94 


0.095 


pump (RSO) -7.53 


16.86 


5.00 


2.94 


0.095 


dump (RSO) 1427.61 


98.32 


1.00 


14.72 


0.474 
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Time (ps) 

FIG. 12: Schematic diagram of the laser pulse sequence con- 
sidered in the present study. The pump and dump pulses are 
centered at 12.0 ps and 25.8 ps, respectively, with a time de- 
lay of 13.8 ps (note that 1 x 10~^ atomic units of electric field 
correspond to 5.14 MY /m). 



D. Photoassociation with a short laser pulse 

The pump-dump scheme considered here is shown 
schematically in Fig. [T^l Each laser pulse is assumed 
to be transform-limited, 



(8) 



with a Gaussian profile f{t) = exp — a(< — ic)^. The pulse 
has central frequency u;l/{2'k) and maximum field am- 
plitude Eq. The full width at half maximum (FWHM) 
of the intensity profile Elf{tf is tl = {2\a2/ a)'^/"^ . 
The corresponding FWHM frequency width is 5u)l{2Tr) = 
2ln2/ (tttl). The Gaussian envelopes are centered at 
times t = tp and id for the pump and dump pulses re- 
spectively. 

The parameters used for the laser pulses in this study 
are given in Table IIIII We choose the parameters of the 
pump pulse to excite a few excited-state vibronic levels 
v' close to the maxima in Fig. [10] where nonadiabatic 
coupling is strongest. The central frequency is chosen to 
be resonant with the level v' = 365, with a laser detuning 
from the atomic line i5|f = 7.53 cm~^ in the RSO case 
and 10.01 cm~^ in the ASO case. The temporal width 
of the pulse is chosen to give an energy spread that will 
excite 5 to 10 levels near v' = 365. The intensity is chosen 
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FIG. 13: The total populations Pi{t) of the excited A and 
b'^n electronic states during the RSO pump pulse, relative to 
the ground-state population. The fractions of the population 
that remain in the excited states after the pulse are shown in 
the inset. 



to obtain maximum excitation of the resonant level. 

The initial state '^{t = 0) is chosen to be a box- 
quantised eigenfunction of the field- free Hamiltonian that 
represents an s-wave scattering state of the ground-state 
potential curve with energy corresponding to a tempera- 
ture T = 16.5 fiK. 

The wavepackct at time t may be written 



^^{R,t)ij}, i = g, 1,3, 



(9) 



where g indicates the ground electronic state and 1 and 
3 refer to the singlet and triplet excited states. The 
wavepackct dynamics are analyzed by studying the evo- 
lution of the population on the respective states. 



P,{t) = 



(10) 



We also project the wavepackct onto the vibronic wave- 
functions of the ground and coupled excited states. 



i=l,3 



(11) 

(12) 



where the vibronic function on the coupled excited states 
is given by Eq. ([7]). Since the amplitude of ^'(t = 0) at 
short range depends on the size of the box used (extent 
of the FGH grid, i?max), the absolute values of the pop- 
ulations and projections decrease approximately linearly 

with i?,„ax- 

The overall populations of the two electronically ex- 
cited states, Pi{t), arc shown in Fig. [13] for the RSO 
case. They reach a maximum near Tp = 12.0 ps. Before 
the end of the pump pulse, most of the population re- 
turns to the initial continuum state. This corresponds to 
coherent transients [s^ which are off-resonant and can 
be excited only during the pulse. Only a small amount 



('--^ 5 X 10"'') remains in the vibronic levels of the cou- 
pled excited states (as shown in the inset of Fig. fTS]) . 
This population continues to oscillate between the two 
electronically excited states due to the nonadiabatic cou- 
pling. Such oscillations have also been observed for Rb2 
[2l|, l53|. The resuhs for the ASO case are qualitatively 
similar, but reduced by about a factor of 2 because of the 
larger detuning from the atomic line. 



Detailed information about the population of the indi- 
vidual vibronic levels w' can be obtained from the projec- 
tion of the excited wavepackct onto excited-state vibronic 
levels PJ' (t). Many vibronic levels are excited transiently 
during the pulse, but only levels with v' = 363 to 370 re- 
main significantly populated after the pulse. The final 
populations of the individual near-resonant vibronic lev- 
els are shown in Fig.[T31 with the time-dependence of the 
populations during the pulse shown in the insets. The 
final populations peak near v' = 365 in both cases, but 
with variations arising from differences in FCFs and for 
levels that are further off-resonance. It is noteworthy 
that v' = 365 itself is almost unpopulated in the ASO 
case: it has a very small FCF because it resides on the 
upper adiabatic curve and has very little probability den- 
sity near the outermost turning point. In the RSO case, 
all the levels have significant upper-state character: in 
this case the anomalously low population of v' = 367 oc- 
curs because the maximum of the last peak of the eigen- 
function corresponds to a node in the initial ground-state 
wavefunction at i? 34.7 oq. 



In order to determine how many molecules are formed 
per photoassociation pulse, it is necessary to average over 
the final excited-state populations obtained from all ther- 
mally populated initial scattering states [131 ■ The aver- 
aging procedure relates the scattering states employed in 
the calculations to the actual volume of the trap. The 
absolute number of molecules is limited by the probabil- 
ity density of atoms pairs near the Condon radius, i.e. 
by the population within the 'photoassociation window'. 
While the exact value depends on the details of the po- 
tential, and in particular on the scattering length, an 
order-of-magnitude estimate can be obtained by compar- 
ing simulations for RbCs and Rb2 using a similar scat- 
tering length, the same size of the grid and an identical 
initial state. In order to compare to the results for Rb2 
, we have repeated the simulation for RbCs shown 
in Fig. [13] using a grid of 20000 ao and an initial state 
with scattering energy corresponding to 100 fiK. The fi- 
nal cxcited-statc norm of 1.1 x 10"'' is about two orders 
of magnitude smaller than for Rb2 (cf. Table I of Ref. 
^ih . where about one molecule per pulse is predicted for 
10^ atoms at a density of 10"^° cm^^ The reduced 
photoassociation yield of a potential compared to 

a 1/i?'^ potential requires a significantly higher density 
and/or a larger number of atoms to produce a compara- 
ble population of ground-state molecules. 
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FIG. 14: The final projections (t) of the wavepacket onto the excited-state vibronic eigenfunctions as a function of their 
binding energies. Significant populations of molecules are formed in vibronic levels v' — 363 to 370. The time-dependence 
during the pulse is shown in the insets for the levels with the largest final population. Left-hand panel: results obtained with 
asymptotic spin-orbit coupling (ASO) Right-hand panel: results obtained with ii-dependent spin-orbit coupling (RSO). 



E. Formation of ultracold ground-state RbCs 
molecules 

Some very loosely bound ground-state vibrational lev- 
els are populated directly by the pump pulse in processes 
involving more than one photon, but a dump pulse is re- 
quired to populate more deeply bound levels. The dump 
pulse is activated after an appropriate time delay, when 
a substantial amount of the wavepacket has reached the 
short-range region where there is significant overlap with 
deeply bound levels of the ground electronic state. Since 
coherent effects between the pump and dump pulses can 
be neglected, the pump and dump pulses are treated sep- 
arately. The excited-state population is normalized to 
one after the pump pulse in order to have a direct cor- 
respondence between the populations calculated and the 
probability of formation of ground state molecules. 

It is useful to inspect the probability density distri- 
bution of the excited-state wavepackets as a function 
of time. Snapshots of the A^S"*" component of the 
wavepackets for a few representative times are shown in 
Fig-HSfor both the ASO and RSO cases. The wavepacket 
is initially quite similar in the two cases and starts to 
move inwards under the influence of the long-range po- 
tential, slightly faster in the ASO case because of the 
larger detuning from the atomic line. However, a major 
difference arises after the wavepackets reach the crossing 
point between the diabatic curves. The ASO wavepacket 
remains almost entirely on the lower adiabatic curve, 
and thus has mostly triplet character inside 10 ao. By 
contrast, the RSO wavepacket is made up of levels of 
mixed upper and lower-state character and there is a 
large buildup of singlet probability density near the in- 
ner turning point. This reaches a maximum at about 
26 ps. Movie of the propagation of the A^I]+ and b-^II 
components of the wavepackets are. 

The parameters of the dump pulse are optimized by 



considering the overlap of the excited-state wavepacket 
with the vibrational levels of the ground electronic state, 
Pg (t). Fig. 1181 shows this overlap function of 

time and binding energy. It may be seen that the RSO 
wavepacket has significant overlap with vibrational lev- 
els of the ground electronic state bound by up to about 
1500 cm~^. The overlap is about 2 orders of magnitude 
smaller in the ASO case because of the lack of popula- 
tion near the inner turning point. This arises from the 
weaker nonadiabatic coupling in the ASO case, which 
in turn arises from the stronger spin-orbit coupling near 
the crossing point. The time-dependence of the overlap 
is also rather different in the ASO case, both because 
of the larger detuning and because any population that 
does arise on the upper adiabatic curve is trapped there. 

After examining the overlaps for different ground- 
state vibrational levels, we chose to de-excite the RSO 
wavepacket to level v" = 54 with binding energy 1435.13 
cm~^ (shown in Fig. [1]). A narrow-bandwidth pulse with 
temporal width 1 ps, detuned 1427.61 cm~^ to the blue 
of the pump pulse, is employed to achieve transfer into 
this single vibrational level (because the excited-state 
wavepacket itself has a binding energy of 7.5 cm~^. As 
may be seen in Fig. \W\ the overlap function for this par- 
ticular vibrational level has a local maximum at t ^ 25.8 
ps. The time delay between the pump and dump pulses is 
therefore chosen to give the nraximum of the dump pulse 
at id = 25.8 ps, corresponding to a time delay of 13.8 ps, 
as shown in Fig. 1121 This is in reasonably good agree- 
ment with the classical prediction of 12 ps for a detuning 
of 7 cm- 1 in Table H 

The final ground-state population Pg'^{t) is shown in 
Fig. [201 as a function of the dump pulse energy. The 
renormalisation applied after the pump pulse removes 
any population that is already in the ground electronic 
state before the dump pulse. Since the dump pulse is res- 
onant only with the bound v" = 54 ground-state level. 
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FIG. 15: Snapshots of the A^E"*" component of the wavepackets at different times after the pump pulse. Left-hand panels: 
results obtained with asymptotic spin-orbit coupling (ASO). Right-hand panels: results obtained with _R-dependent spin-orbit 
coupling (RSO). 



FIG. 16: Movie showing the A^E"^ and b'^II component of 
the wavepackets as a function of time for the asymptotic spin- 
orbit coupling (ASO) coupling function. [Animated gif file 
movieASO.gif] 



FIG. 17: Movie showing the A^E''' and b"^!! component of 
the wavepackets as a function of time for the i?-dependent 
spin-orbit coupling (RSO) coupling function. [Animated gif 
file movieRSO.gif] 



and the ground-state vibrational spacing of about 37 
cm~^ is much larger than the bandwidth of the pulse, 
the population of the v" = 54 level is the same as the to- 
tal ground-state population. The maximum population 
transfer is obtained for a pulse of energy ~ 9 ^J. We 
repeated the dump calculation for the ASO wavepacket 
(with a modified laser frequency and a slightly differ- 
ent time delay to correspond with the maximum in Fig. 
fTS]) and verified that the population of ground-state 
molecules is about two orders of magnitude smaller. 

Ground-state RbCs molecules bound by as much as 
1500 cm~^ can be formed because there is good Franck- 
Condon overlap of the excited-state wavepacket with 
deeply bound vibrational levels. Similar studies on the 
homonuclear Rb2 molecule [U have shown that ex- 
cited Rb2 molecules can be efficiently dumped to create 
ground-state vibrational levels bound by only about 10 
cm~^. For Rb2 the wavepacket showed significant peak- 
ing of singlet population at the outer tuning point of the 
upper adiabatic curve but not at its inner turning point 
as in RbCs. The time delay used for Rb2 was 81.5 ps 
with a detuning of 4.1 cm~^ from the atomic line. 

As commented above, the levels of Rb2 do show strong 
nonadiabatic mixing between the upper and lower adia- 
batic states. In view of this, it is important to understand 



why RbCs shows much more buildup of singlet charac- 
ter near the inner turning point than Rb2. The reason 
for this is that Rb2 has an R^"^ potential at long range, 
so that the vibrational wavefunctions of states near dis- 
sociation are more strongly concentrated at the outer 
turning point. Because of this, even population that is 
transferred to the upper adiabatic state by nonadiabatic 
coupling has a smaller amplitude near the inner turn- 
ing point and is less effective in allowing deexcitation to 
deeply-bound states. 

We also note that, because of the small FC factors 
for deexcitation, only about 3% of the excited-state 
RbCs population can be transferred to the desired deeply 
bound ground-state vibrational level, whereas almost 
50% of the excited Rb2 population can be transferred 
to the ground state using much less dump pulse energy 
(Fig. 11 of 2l|)- The overall efficiency of the scheme 



considered here is limited by the fraction of the initial 
ground-state population that is transferred to the ex- 
cited state. One possible way to enhance it would be 
by Feshbach-optimised photoassociation [ssj . 



VI. CONCLUSIONS 

We have studied the possibility of using photoassoci- 
ation with laser pulses to produce deeply bound RbCs 
molecules in the electronic ground state. We have also 
explored the differences between photoassociation pro- 
cesses for RbCs and Rb2 in order to understand the con- 
sequences of the different long-range potentials in excited 
states of heteronuclear and homonuclear molecules. 

The difference between potentials for heteronu- 
clear molecules and potentials for homonuclear 
molecules produces several important effects. First, the 
Franck-Condon factors for excitation are smaller in the 
heteronuclear case, though they die off more slowly with 
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FIG. 18: The time-dependence of the projections Pg (t) of the wavepackets onto the vibrational levels of the ground electronic 
state. Left-hand panel: results obtained with asymptotic spin-orbit coupling (ASO). Right-hand panel: results obtained with 
7?-dependent spin-orbit coupling (RSO). 
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FIG. 19: The projection of the excited-state wavepacket onto 
the v" — 54 vibrational level of the ground electronic state. 
The time indicated by the arrow is the maximum of the dump 
pulse (td). 
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FIG. 20: Fraction of the excited-state population transferred 
to the v" = 54 level of the X ^ ground state as a function 
of the dump pulse energy for a time delay of 13.8 ps. 



detuning from the atomic line. Secondly, the time taken 
for a wavepacket produced on the upper state to evolve 
from long to short range is significantly shorter for a het- 
eronuclear molecule than for a homonuclear molecule of 
similar mass at the same detuning. We give a simple 
semiclassical expression relating the time delay to the 
detuning and the coefficients governing the long-range 
potential. 

We have also explored the dependence of Franck- 
Condon factors on the scattering lengths for both the 
ground state and the excited state. The Franck-Condon 
factors oscillate as a function of detuning, with faster os- 
cillations for RbCs than for Rb2. Adjusting the ground- 
state scattering length alters the positions of the peaks 
in the Franck-Condon factors as a function of detuning, 
while adjusting the excited-state scattering length leaves 



the oscillations unchanged but shifts the vibrational lev- 
els within them. 

We have carried out wavepacket calculations to ex- 
plore the quantum dynamics of RbCs on 0+ excited states 
formed by coupling the A and b ^11 electronic states. 
The diabatic potential curves for the A and b states cross 
near R = 10 oq. The dynamics are strongly affected by 
the magnitude of the spin-orbit coupling. If the spin- 
orbit coupling function is held constant at its asymptotic 
value, the dynamics takes place almost independently on 
the upper and lower adiabatic curves. Under these cir- 
cumstances the wavepacket has almost entirely triplet 
character at short range and the Franck-Condon factors 
for deexcitation to deeply-bound levels of the ground 
state are very poor. However, if the spin-orbit coupling 
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is given a more realistic experimentally derived form |37| 
that has a smaller value near the crossing point, the dy- 
namics are strongly nonadiabatic. A substantial part of 
the population is transferred to the upper adiabatic state, 
which has singlet character at short range. This allows 
efficient deexcitation to levels of the ground electronic 
state bound by up to 1500 cm~^. 

The behavior observed for RbCs may be contrasted 
with that for Rb2 [HI, • R.b2 has only slightly weaker 
nonadiabatic coupling but deexcitation is favored to 
ground-state levels bound by up to only 10 cm^^. The 
difference in this case arises from the different long-range 
potentials: for Rb2 there is an potential at long 
range and the density distribution for the upper-state 
levels is dominated by the outer turning point of the up- 
per adiabatic curve. From the outer turning point, the 
Franck-Condon factors favor deexcitation to relatively 
weakly-bound levels of the ground electronic state. 

An important new insight from the present work is that 
the combination of strong nonadiabatic coupling with a 
potential produces mixed vibrational levels with 
significant singlet density both at long range (which facil- 
itates initial photoassociation) and near the inner turning 
point of the singlet state (which facilitates deexcitation to 
form deeply bound ground-state levels). This allows 'i?- 
transfer' in a way which does not occur for homonuclear 
species. This phenomenon is also expected to be present 
in other heavy heteronuclear dimers such as KRb or KCs. 
It is easily identified spectroscopically in the level spac- 
ings or rotational constants and will play an important 
role that is not confined to photoassociation, for exam- 
ple in stimulated Raman adiabatic pumping (STIRAP) 
experiments to form deeply bound states. 

Our study of pump-dump photoassociation was mo- 
tivated by the search for a coherent scheme for form- 
ing deeply bound molecules where the dump pulse does 
not destroy the molecules that were created by the 
pump pulse. STIRAP is another candidate for this, and 
has recently bee n applie d successfully to transfer near- 
dissociation Cs2 [l6l . Isg} and KRb [l^, molecules to 



deeply bound states. However, STIRAP requires a few 
well-separated levels in order to fulfill the adiabaticity 
condition and this requirement is not met by the quasi- 
continuum of a MOT. This might be circumvented by 
applying STIRAP after a photoassociation pump pulse, 
in order to optimize the dump step. In other words the 
pump pulse prepares a coherent state which in a two- 
level picture replaces the upper level, while the target 
ground-state vibrational level plays the role of the lower 
level. Several realizations of adiabatic passage in two- 
level systems exist, such as Stark-chirped rapid adiabatic 
passage (SCRAP) [H, [5§| and retroreflection-induced 
bichromatic adiabatic passage (RIBAP) [13, [Hj. The 
underlying concept of these schemes is to induce a cross- 
ing of the two levels and an adiabatic or diabatic passage 
of the induced crossing. In our case the adiabatic passage 
would need to proceed on a timcscalc long compared to 
the vibrational motion of the wavepacket (about 40 ps) 
and short compared to spontaneous emission (about 30 
ns). However, to see whether replacing the picosecond 
dump pulse by an adiabatic passage scheme is feasible 
and whether it allows the dump efficiency to be increased 
above the 3% reported here is beyond the scope of the 
present study. 
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